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ABSTRACT 

The Infrared Astronomical Satellite (IRAS) had a tremendous impact on many areas of modern 
astrophysics. In particular it revealed the ubiquity of infrared cirrus that are a spectacular mani- 
festation of the interstellar medium complexity but also an important foreground for observational 
cosmology. With the forthcoming Planck satellite there is a need for all-sky complementary data sets 
with arcminute resolution that can bring informations on specific foreground emissions that contam- 
inate the Cosmic Microwave Background radiation. With its ~ 4' resolution matching perfectly the 
high-frequency bands of Planck, IRAS is a natural data set to study the variations of dust properties 
at all scales. But the latest version of the images delivered by the IRAS team (the ISSA plates) 
suffer from calibration, zero level and striping problems that can preclude its use, especially at 12 
and 25 /im. In this paper we present how we proceeded to solve each of these problems and enhance 
significantly the general quality of the ISSA plates in the four bands (12, 25, 60 and 100 /im). This 
new generation of IRAS images, called IRIS, benefits from a better zodiacal light subtraction, from 
a calibration and zero level compatible with DIRBE, and from a better destriping. At 100 /im the 
IRIS product is also a significant improvement from the Schlegel et al. (1998) maps. IRIS keeps the 
full ISSA resolution, it includes well calibrated point sources and the diffuse emission calibration at 
scales smaller than 1° was corrected for the variation of the IRAS detector responsivity with scale and 
brightness. The uncertainty on the IRIS calibration and zero level are dominated by the uncertainty 
on the DIRBE calibration and on the accuracy of the zodiacal light model. 
Subject headings: ... 



1. INTRODUCTION 

From January to November 1983 the Infrared As- 
tronomical Satellite (IRA S), a joint project of th e US, 
UK and the Netherlands f Neugebauer et al.l ll984'l. per- 
formed a survey of 98% of the sky at four wavelengths 
: 12, 25, 60 and 100 /im. IRAS lead to numerous scien- 
tific discoveries spanning a broad range of astrophysical 
subjects, from comets, circumstellar disk to interacting 
galaxies. The satellite was designed to optimize the re- 
liablility of point sources detection and photometry; one 
of the great legacies of IRAS is certainly its catalog of 
more than 250 000 point sources. On the other hand, the 
relative stability of its detectors also allowed mapping of 
extended emission. In fact IRAS made a significant con- 
tribution to our understanding of Galactic diffuse emis- 
sion by revealing the interstellar dust emission of infrared 
cirrus, that can b e observed in any direction on the sky 
lILow et al.lll98l . 

A first set of extended emission maps (the SkyFlux at- 
las) was released in 1984 and 1986 along with the IRAS 
Point Source Catalog. It was then clear that significant 
improvements of the sensitivity and photometric accu- 
racy could be obtained based on the acquired knowledge 
of the instrument. The IRAS data were later reprocessed 
based on the knowledge of the instrumental response 
available at that time but also by applying a zodiacal 
light correction and a destriping method. This second 
generation processing, which increased the sensitivity by 
a factor up to five, lead to the the IRAS Sky Survey At- 



las (ISSA) published in 1991 and 1992. The IRAS Sky 
Surve y Atlas Explanatory Supplement (Wheclock et al] 
Il993|) describes in detail the processing and analysis of 
this atlas. 

Since its publication, the ISSA has been a tremen- 
dously useful data product, used to study various aspects 
of Galactic (e.g. Boulanger & Perault (1988)) and extra- 
galactic (e.g. Miville- Deschenes et al.l (|2002b|)) diffuse 
emission. The ISSA plates have become an essential data 
set for any multi-wavelength analysis of the interstellar 
medium. Unfortunately, even after the heroic processing 
done by the IRAS team, the ISSA plates still suffer from 
defects (striping, calibration, zero level, zodiacal light) 
that can limit significantly their use. In the context of the 
numerous ongoing and future infrared/submm/mm mis- 
sions (Spitzer, Herschel, Planck, JWST, ASTRO-F...), 
IRAS data could still be highly valuable, even with their 
relatively coarse angular resolution (^^4'). This was our 
main motivation to reprocess the ISSA plates of all IRAS 
bands, in order to improve their sensitivity and absolute 
calibration. 

This new generation of IRAS maps, a product named 
IRIS (for Improved Reprocessing of the IRAS Survey), 
will be of particular interest for the Planck mission ded- 
icated to the study of the Cosmic Microwave Back- 
ground (CMB). Planck will perform an all-sky survey 
in the 350 /im-10 mm wavelength range with an an- 
gular resolution in the high frequency bands compara- 
ble to IRAS. One of the main anticipated challenges 
of this mission will be the separation of all the com- 
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ponents, including dust, that contribute to the sig- 
nal in this wavelen gth range (as it is for WMAP, 
ijBennett et alJ l2003|) '). Several studies of dust emis- 
sion at high Galactic latitudes showed large variations of 
dust properties at all scales, even in the diffuse interstel- 
lar medium (Boul^ne rer et al. 1990: Laeache et al.. ,199&: 
Bernard et all 119991: iMiville-Deschenes et alJ l2002al: 
Stepnik et 311120031). T o imp rove o n the dust mod- 
els available (e.g. iFink beiner et alJ ljl999|) ') and bring 
them to an accuracy level required for Planck, dust 
evolution should be understood in a more general way, 
including the properties of the smaller grain popula- 
tions that give some informations on the evolutionary 
stage of dust (coagulation, fragmentation) and on the 
properties of the bigger grains that contribute to the 
submm/m m emission. Smal l grains have also been 
suggested l|Draine fc LazarianI [l998) to be re sponsible 
for the anomalous microwave emission (e.g. iLaeachel 
((2003)). In this context, well calibrated, arcminute reso- 
lution and all-sky dust emission maps at 12, 25, 60 and 
100 fim would be of tremendous help to understand and 
model dust emission in the Planck bands. This is one of 
the main motivations of the present work. 

The present paper describes in detail the procedure 
we followed to reduce significantly the striping of the 
ISSA maps and to rescale them to have an absolute 
calibration coherent with both the DIRBE and IRAS 
Point Source Catalog photometry. In § |21 we describe 
in detail the ISSA product. Our reprocessing of the 
ISSA plates, which include an automatic deglitching, 
an improved destriper and an absolute calibration and 
zero level correction, is described in § 13 and § 0| 
The overall properties of the IRIS product are dis- 
cussed in § |S1 The whole IRIS data set is available at 
|http://www.cita.utoronto.ca~mamd/IRIS and at 
[http: / /w ww. ias.u-psud.fr / iris / 

2. THE IRAS SKY SURVEY ATLAS 

2.1. Processing of the ISSA plates 

During its 10 months operation period, IRAS made 
two surveys of 98% of the sky and a third one of 75% of 
the sky. To confirm point source detection, each survey, 
called an HCON for Hours CONfirmation, consisted of 
two coverages of the sky separated by up to 36 hours. 
The first two HCONs (HCON-1 and HCON-2) were done 
concurrently with the second one behind the first one by 
a few weeks. The third survey (IICON-3) began after the 
first two were completed. Due to exhaustion of the liquid 
helium supply the third HCON could not be completed 
and it thus covers only 75% of the sky. 

The ISSA is a set of 430 fields of 500 x 500 pixels of 
1.5'. Each field is a 12.5° x 12.5° region centered every 
10° in declination. All HCONs represent the sky sur- 
face brightness at the 4 wavelengths from which an esti- 
mate of the zodiacal light was subtracted. The HCONs 
were destripped, checked for zero level stability and visu- 
ally examined for artifacts and glitches. The processed 
HCONs were then co-added, using sky coverage maps, 
to produce one final map (named HCON-0) per field and 
per wavelength. 

2.2. Diffuse emission photometry 

Special care has been taken to have a consistent diffuse 
emission calibration through the Atlas. But, as the IRAS 



mission was designed to provide absolute photometry for 
point sources only, the ISSA images give only relative 
photometry on spatial scales larger than ~5' and can- 
not be used directly to determine the absolute surface 
brigthness of diffuse emission. The uncertainty on the 
zero level is dominated by uncertainties in the zodiacal 
light model at 12 and 25 fim, and by detector drifts at 
60 and 100 /im. 

As it was designed to study the cosmic infrared back- 
ground, the DIRBE (Diffuse Infrared Background Exper- 
iment) experiment on board the COBE (Cosmic Back- 
ground Explorer) satellite has made absolute surface 
brightness measurements of the sky at scales larger than 
0.7°. Moreover the DIRBE data were correct ed with a 
better zodiacal light model (Kclsall et al. 1998) than the 
one used for the ISSA. Fortunately DIRBE and IRAS 
had very similar filters and therefore the DIRBE data 
can be used to determine the absolute IRAS calibration 
at scales larger than 0.7°. Some efforts have been put 
in the comparison between DIRBE data and individual 
HCON images. Wheclock et al. (1993) showed that the 
amplitude of the 60 and 100 /im intensity fluctuations are 
overestimated in the IS SA maps with respect to DIRBE. 
iWheelock et all ()1993(1 suggest a linear relation between 
the IRAS and DIRBE surface brightness that also takes 
into account the difference in zero level: 

^^^^(a, <5) = Ga X ^-^(a, S) + B^ia, 6). (1) 

The gain correction factors (Ga) given in Tabled were 
originally obtained on spatial scales of the order of an 
IRAS scan 8°). We have performed a detailed com- 
parison of the DIRBE and IRAS data on numerous re- 
gions covering a broad range of brightness and we con- 
clude that these factors vary with angular scale and 
brightness. This issue is discussed in detail in i) l4.3| Re- 
gardin g the zero point of the ISSA maps, iSchlegel et alJ 
1)1998(1 showed that the 100 /xm offset Bioo{a,d) varies 
significantly across the sky. Our analysis show that this 
is the case at all wavelengths (see § 14.5)1 . 

2.3. Point Source Photometry 

By st aring at point sources , the IRAS team clearly 
showed f Wheelock et al.l 11993)) that the instantaneous 
(AC) and long-term (DC) responsivity of the IRAS de- 
tectors were different. For a survey mission like IRAS, 
this introduces a responsivity difference between point- 
like object and diffuse emission. The ISSA plates are 
DC calibrated which means that point source fluxes es- 
timated on a ISSA plate have to be corrected for the 
DC/ AC responsivity difference (see Table Even if 
the ISSA product was not optimized for point sources 
fWhcclock et al.l 1)19931 ) showed that, when corrected for 
the appropriate DC/ AC factor, point source fluxes mea- 
sured on ISSA plates are photometrically consistent with 
the IRAS Point Source Catalog to within 10%. We have 
verified this by comparing the DIRBE and ISSA (con- 
volved by the DIRBE beam) fluxes of several stars. We 
also found a general good agreement between these two 
data sets. 

2.4. Post-processing of the ISSA plates 

Beside the responsivity variation and zero level prob- 
lems, the ISSA plates have also relatively strong stripes. 
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TABLE 1 

Compilation of IRIS related quantities 



Wavelength (/im) 



Quantity 


12 


25 


60 


100 


1. 


IRAS resolution (arcmin) 


0.75 X 4.5 


0.75 X 4.6 


1.5 X 4.7 


3.0 X 5.0 


2. 


IRIS resolution (arcmin) 


3.8 ±0.2 


3.8 ±0.2 


4.0 ±0.2 


4.3 ±0.2 


3. 


ISSA noise level (MJy/sr) 


0.04 ±0.01 


0.05 ±0.02 


0.04 ±0.01 


0.07 ±0.03 


4. 


IRIS noise level (MJy/sr) 


0.04 ±0.01 


0.05 ± 0.02 


0.03 ±0.01 


0.06 ±0.02 


5. 


DC/AC 


0.78 


0.82 


0.92 


1.00 


6. 


Gx 


1.06 ± 0.02 


1.01 ± 0.02 


0.87 ±0.05 


0.72 ±0.07 


7. 


DIRBE S(gain) (%) 


5.1 


15.1 


10.4 


13.5 


8. 


DIRBE S(offset) (MJy/sr) 6.0 X 10^5 


8.3 X 10-^ 


0.027 


0.027 



Note. - Row 1: Effec tive instrumental resolution of the IRAS detectors IWheelock et al1ll993|1 . Row 2; Angular resolution of the ISSA 
and IRIS plates (see S 15.21 . Row 3-4: The ISSA and IRIS noise level given here are representat ive o f the median instrumental noise of the 
whole plates. The noise level varies f rom one plate to the o ther as a function of coverage (see i; l5.1i . Row 5-6: D C/AC and IR AS/DIRBE 
corrections given by the IRAS team IWheelock et al1ll993l) . Row 7-8: DIRBE gain and offset uncertainties from lHauser et all 11998) . 



related to the scanning on the sky, that participate sig- 
nificantly to the noise. The main goals of the present 
work were first to reduce the amplitude of those stripes 
to lower the noise level and second to correct the ISSA 
plates calibration to match both the DIRBE and the 
Point Source Catalog photometry. 

In the next sections we describe how we have addressed 
all of these problems. First we discuss how we have 
dcglitched automatically the individual HCONs fS 13. 
Then we describe our destriping algorithm I3.2|l . Fi- 
nally we describe how we use the comparison of the IRAS 
and DIRBE data to establish a calibration and a zero 
level correction that is coherent both for point sources 
and diffuse emission (§0)1. 

Through out the processing we have taken great care 
to work on individual HCONs to improve the destriping 
of the maps and to be able to estimate the noise level in 
the final averaged maps. This specific point is discussed 
in §0 

3. CLEANING THE ISSA PLATES 
3.1. Deglitching 

Before co-adding the individual HCONs, the IRAS 
team inspected them visually to identify artifacts like 
moving objects (e.g. asteroids) or residual glitches. To 
allow for further validity checks of point sources seen 
(or not) in the co-added maps, the IRAS team decided 
to provide the individual HCONs that include those ar- 
tifacts. As we are working directly on the individual 
HCONs we had to redo the deglitching but we took ad- 
vantage of the fact that the co-added maps (HCON-0) 
have already been deglitched. 

First, using the method described in § ^ we identi- 
fied point sources both in the individual HCONs and in 
the co-added map. In a given map the identified pixels 
were then gathered together into regions by consecutively 
labelling all the detecting regions with a unique region 
index ^. Then we computed the central position (taken 
as the position of the brightest pixel) and the maximum 
radius of each region. For each region found in the in- 
dividual HCON (1,2 or 3) we derived the distance to all 
regions found in the co-added map. Each HCON region 
with no counterpart in the co-added map at distance 

^ We use the label_region function of IDL 



lower than the sum of the two regions maximum radius 
is considered as a glitch. The corresponding pixels in 
the HCON are then set to the background level at that 
position. 

This automatic method worked very well to identify 
automatically residual glitches and moving objects on 
most of the ISSA plates. It was not as performant in 
very bright regions were the diffuse emission varies sig- 
nificantly at small scales. Therefore we did not deglitch 
regions brighter that 200 MJy/sr at 100 /im. We also 
dealt manually with Saturn due to its large extend in 
the affected HCONs. 

3.2. Destriping 

The ISSA plates have relatively strong stripes that are 
related to the scanning on the sky. We describe in this 
part our destriping algorithm that reduces significantly 
the amplitude of those stripes. The destriping is done in 
Fourier space for each HCON. For typical HCONs, the 
stripes have one or two directions but there exists more 
complicated regions where there are several directions 
(up to 5 for example near the Large Magellanic Cloud). 
In the Fourier space, such stripes give extra-power in 
well-defined azimuthal directions that can be identified 
and corrected (Fig. 

To be able to identify clearly the striping pattern in 
Fourier space, two effects have to be dealt with: point- 
like sources and large scale gradients. Point sources in- 
troduce contaminating power in Fourier space in every 
direction and with a flat power spectrum. Their main 
effect is to lower significantly the contrast of stripes in 
Fourier space, especially at large wavenumbers k. On the 
contrary, a large scale gradient across an image produces 
extra-power in privileged directions in Fourier space (hor- 
izontal and vertical) that mimics stripes. To get rid of 
this effect it is quite common to apodize the image (see 
ISchlcgel et al. (1998) for example) but this destroies the 
information on the edge of the image. 

To get round these difficulties we used a wavelet de- 
composition ("a trou" algorithm) of every HCON and 
worked only on an image that is the sum of the first five 
scales (2, 4, 8, 16 and 32 pixels of 1.5'). This filtered map 
contains only small scale structures, usually dominated 
by the stripes, and is thus free of any large scale gradient 
(Fig. nj. To increase the contrast of the striping pattern 
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we removed point sources at each scale of the wavelet 
decomposition prior to summing them; we cliped all the 
wavelet coefficients higher than 7a, where a is computed 
at each scale. The original HCON was then decomposed 
into three components: a large scale emission map (scales 
larger than 32 pixels), small scale emission map (scales 
smaller than 32 pixels) and a point sources map. The 
following destriping algorithm was only applied on the 
small scale emission map. 

Each small scale emission map was Fast Fourier Trans- 
formed (FFT) and the averaged magnitude of the FFT 
was computed as a function of angle 9 in 360 bins of 1° 
wide. We identified the bins with averaged magnitude 
greater than Ntr and replaced them by the averaged az- 
imuthal magnitude. This is equivalent to replace the 
identified bins by the power spectrum of the clean bins. 
This procedure was repeated twice on artificially rotated 
FFT pattern (to detect stripes with direction falling be- 
tween two bins). During this process, we progressively 
decreased N from 7 to 5. We then made the inverse 
transform of the cleaned FFT and added the result to 
the large scale emission and point source maps. 

This procedure can only be applied on fully defined 
HCONs. For each plate, the missing data were itera- 
tively replaced by the destriped data of the companion 
HCONs. When large regions have not been observed by 
the 3 HCONs, we used the DIRBE data, corrected for 
the appropriate gain and offset values. 

Our procedure has been extensively tested on sim- 
ulations. We simulated typical interstellar dif- 
fuse regions using fractional B rownian motion images 
ijMiville-Deschenes et al.l 120031) with a similar power 
spectrum to what is observed in the ISM (spectral in- 
dex of -3 (|Gautier et al.lll992D 'l. We added to this map 
typical random noise together with stripes that were sim- 
ulated in Fourier Space (note that we cannot inject in 
our simulations empirically derived stripes from real data 
since such maps, obtained by wavelet filtering real data 
visually dominated by stripes, are also contaminated by 
the instrumental noise). The simulated map were then 
processed through the whole destriping procedure. An 
example of such a simulated map for a typical high lati- 
tude region at 12 /xm (where the amplitude of the stripes 
are of the same order than the signal) is shown on Fig. 
El together with the result of the destriping algorithm. 
We computed the power spectrum of the input noisy in- 
terstellar map and that of the recovered map. In such 
regions, we modified the power spectrum at small scales 
by only about 3%. This demonstrates that we have su- 
pressed the contribution of the stripes (that was four 
times higher than the signal on the power spectrum at 
small scales in this case) without adding extra-noise. In 
brighter regions, the contribution of the stripes becomes 
smaller compared to the galactic emission, leading to a 
still better recoverage of the simulated map. These sim- 
ulations clearly show the robustness and efficiency of our 
destriping process. 

4. CALIBRATING THE ISSA PLATES 

4.1. Infrared photo- detectors 

iWheelock et al.l l|1993|) showed that the responsivity of 
all IRAS detectors are affected by memory effects that 
introduce a calibration difference between diffuse emis- 
sion and point sources. Since IRAS, and especially with 
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Fig. 1. — Illustration of the destriping algorithm. The image 
(plate number 83 at 60 fira) is choosen such that it contains point 
sources, diffuse emission and undefined regions. Top-left: HCON; 
Top-right: wavelet filtered HCON; Middle-left: FFT of the fil- 
tered HCON; Middle-right: modified FFT; Bottom-left: destriped 
HCON; Bottom-right: ratio of the power spectrum of the HCON 
and the destriped HCON. 



the ISO mission, several studies have been done to un- 
derstand the response of infrared photo conductors like 
the IRAS 60 and 100 /xm Ge:Ga de tectors. Using solid 
state and semi conduct or phy s ics. [ Vinokurov fc Foukd 
1)19911) and |^ks & Sc hubertI l)1995|) have studied the 
response of such detectors. They showed that the in- 
stantaneous response depends on the illumination and 
on the history of the detector. A simplify version of 
their m odel was succ essfully applied to ISOCAM obser- 
vations l|Coulais fc A bergel 20001. This non-linear model 
is rather complex and it is not the purpose of the present 
paper to describe it in detail. But to feel how the re- 
sponse of infrared photo detectors behaves, lets consider 
the simple case of a detector stabilized at flux Jq that is 
suddenly illuminated at flux Ji at time i = 0. In this 
case, the response J{t) of the detector as a function of 
time is described by: 

Jq + [Ji - Jo) exp[-tJi/X) 

The (3 and A parameters depends on the physical proper- 
ties of the detector. An example of this model is shown 
in Fig. where we illustrate what would be the response 
of a detector when the flux is double (Ji = 2 x Jq), for 
three values of the initial flux Jq. We see that the re- 
sponse of the detector is very sensitive to the flux value; 
for lower flux value, the detector takes more time to stabi- 
lize. Therefore, for a given integration time, an intensity 
fluctuation on a low surface brigthness region will be sys- 
tematically underestimated in comparison to a brighter 
one. 

This has an important consequence on survey observa- 
tions like the IRAS mission that scanned the sky with a 
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Fig. 2. — Simulated high latitude cirrus region with noise and stripes. Top-left: cirrus (P{k)oc k~^); Top-middle: random noise; Top- 
right: stripes; Bottom-left: total simulated map; Bottom-middle: destriped map; Bottom-right: power spectrum of the cirrus (continuous 
line), cirrus+noise (dashed line), cirrus-|-noise-|-stripes (diamonds), and the recovered map after the destriping process (dashed-dotted line), 
in arbitrary units. The dashed and dashed-dotted lines are superimposed (the difference is at the order of 3% at very small scales). 



fixed integration time. It is likely that, at low flux value, 
the integration time was shorter than the stabilization 
time of the detector. Therefore the amplitude of the fluc- 
tuations on the sky are systematically underestimated in 
faint regions. This type of memory effect should leave 
another signature on maps with a fixed integration time. 
As the satellite scans the sky it observes structures with 
various angular scales. These structures can be seen as 
flux steps like in Fig. 13 Then, when the satellite goes 
through a small scale structure, it is likely that the ob- 
served flux will never reach the stabilization level. On the 
other hand, for very large scale structures, the detector 
has enough time to stabilize. Therefore, the amplitude 
of the intensity fluctuations are underestimated at low 
flux but also at small scales. 

4.2. Responsivity of the 12 and 25 /jm detectors 

Similar effects were indeed observed by the IRAS 
team. By staring at se veral stars of different brightness, 
IWheelock et al.l 1)199,'^ showed that the responsivity R 
of the 12 and 25 micron detectors depends on time and 
therefore on angular scale, as the satellite is scanning 
the sky at a fixed rate (3.85 arcmin/sec). The responsiv- 
ity was observed to be constant for angular scales larger 
than 2° but it decreased at smaller scales, in accordance 
with the model described earlier. On the other hand 
the IRAS team reports no variations of the responsiv- 
ity as a function of brightness for these detectors. In 
Table |21 w e give the respon sivity variation with angular 
scale from IWheelock et al.l |1993') for the 12 and 25 /im 
detectors but one must be aware that this variation of 
the responsivity with scale was not applied to the ISSA 
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Fig. 3. — Illustration of the response of an infrared photo de- 
tector based on Eg. 1^ I Fouks & Schubert 1993). The response to 
upward and downward flux stops are shown for three values of the 
flux Ji = 2, 6 and 10. To illustrate the behavior of the detec- 
tor when the flux value is double, the lower flux value is Ji/2 for 
each case. The response of the detector is very sensitive to the 
flux value; for lower flux value, the detector takes more time to 
stabilize. 



plates. Instead a constant responsivity factor was used to 
match the responsivity at large scale (i.e. the ISSA plates 
are DC calibrated) which means that the brightness fiuc- 
tuations at scales smaller than ~ 2° are systematically 
underestimated at 12 and 25 /im. 
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TABLE 2 

Responsivity factors 



I 


12 /zm 


25 


I 




60 /im 






100 ^lm 




(degree) 


Ro 


i?.o 


(degree) 


Ro 


A 


B 


Ro 


A 


B 


0.067 


0.78 


0.82 


0.067 


0.92 


0.0 


0.0 


1.0 


0.0 


0.0 


0.1 


0.90 


0.90 


1.25 


1.00 


0.08 


-5.0e-5 


0.96 


0.008 


-4.8e-5 


0.2 


0.92 


0.93 


2.5 


1.01 


0.01 


-4.0e-5 


0.99 


0.011 


-6.5e-5 


0.5 


0.98 


0.97 


5.0 


1.05 


0.014 


-6.0e-5 


1.08 


0.012 


-8.0e-5 


1.0 


1.00 


0.98 


10.0 


1.13 


0.014 


-9.0e-5 


1.25 


0.08 


-5.9e-5 


1.5 


1.00 


1.00 
















2.0 


1.00 


1.00 

















The responsivity R at a given wavelength, scale and brightness is given by Eq. |3] The responsivity at 12 and 25 fim does not depend on 
brightness {A = 0, B = 0). The value of the responsivity at the smallest scale (0.067 degree) is the DC/AC factor (see Table Q needed to 
recover the correct point source calibration. 



4.3. Responsivity of the 60 and 100 /im detectors 

At 60 and 100 //m the situation is complicated by the 
fact that the responsivity depends both on time (i.e. an- 
gular scale) and brightness. Wheelock et al. (1993) did 
not give any prescription to address this problem. To 
evaluate the variation of the responsivity with scale and 
brightness of the 60 and 100 /xm detectors we compared 
the IRAS and DIRBE data using a wavelet analysis. 

To study the variations of the 60 and 100 /im respon- 
sivity we have looked at variations of AIRAS/ADIRBE 
as a function of scale and brightness. For a given diffuse 
structure on the sky the quantity AIRAS/ADIRBE 
represents the ratio of the brightness fluctuation (above 
the background) of that structure seen in the IRAS and 
DIRBE data. If we take the DIRBE calibration as an 
absolute reference, the quantity AIRAS/ADIRBE is 
equal to the IRAS responsivity i? to a brightness fluc- 
tuation. To compute AIRAS/ADIRBE as a function 
of scale and brightness we flrst put the two dataset 
on a common grid and at the same angular resolution. 
For each ISSA plate we have constructed a larger map 
(600 X 600 pixels) by mosaicking neighbouring plates. 
This step is necessary in order to avoid edge effects in 
further convolutions. Next this bigger ISSA plate is con- 
volved by the DIRBE beam. The DIRBE data are then 
projected on the ISSA grid. Finally both maps are con- 
volved by a Gaussian be am of FHWM 40' to smooth out 
the DIRBE pixelisation l|Schlegel et alJl998|) and the ex- 
tra 50 pixels wide edges are removed. At that point the 
angular resolution and size of both maps are of 1 degree 
and 500 x 500 respectively. 

To compute the variation of AIRAS/ADIRBE with 
scale and brightness we used a wavelet transform (the 
"a trou" algorithm - IStarck fc Murtaehl l|1998j) ) that al- 
lows us to study the brightness variations as a function 
of scale. An exemple of such a decomposition is shown 
in Fig. El on a fractional Browni an motion (fBm) im- 
age ifMiville-Deschenes et al.l[2003j) with a similar power 
spectrum to what is observed in the ISM (spectral in- 
dex of -3) and convolved the same way the IRAS and 
DIRBE data were. On the top row of Fig. 0] we show 
two fBm images, mi and m2, where m2 was constructed 
by simply applying a constant multiplicative factor to mi 
(m2 = Q X mi). On the four bottom rows we show the 
wavelet coefficients at four differents scales (1.25, 2.5, 5 
and 10 degrees). The important result here is that the 
slope of the linear regression of the wavelet coefficients 



gives exactly the multiplicative factor Q at each scale. 

An exemple of the same wavelet decomposition for a 
typical 60 /im field is shown in Fig. One would notice 
that the wavelet coefficients of the IRAS and DIRBE 
are very well fitted by a linear fit with relatively small 
dispersion. One noteworthy result here is the increase 
with scale of the slope (AIRAS/ADIRBE) of the linear 
regression between DIRBE and IRAS data. The trend 
seen in Fig.Elis totally representative of what is generally 
observed. 

We have perform this analysis at 60 and 100 /im, on all 
the ISSA plates with no undefined values, to also study 
the variation of R = AIRAS/ADIRBE with brightness 
< I\ >^. The results of this analysis are compiled in 
Figures |S1 and [3 where we show the variation of R with 
brightness at the four scales studied. ^From these two 
figures we found that the IRAS responsivity R increases 
with brightness for the two bands. The increase of R with 
scale described previously is also seen in these two fig- 
ures. In addition the values of R are systematically above 
1.0 at large scale which confirms that IRAS brightness 
fiuctuations are overestimated at large scales in the ISSA 
plates. The values obtained here at scales of 10° and for 
brightness of ^ 3 and ^ 10 MJy/sr at 60 and 100 fim 
are com patible with the IRAS-D IRBE correction factors 
given bv lWheelock et alJ l(T993|) (see Table [TJ. 

The increase of R with scale and brightness indicates 
that the IRAS responsivity is generally lower for small 
scale and faint brightness fluctuations than for large scale 
and large brightness fiuctuations. This behavior is in 
complete accordance with what is expected for typical 
photo-conductors (see ij l4.1|l . 

The dotted lines in Figures and [7| represent the slope 
of the finear regression of IRAS vs DIRBE and DIRBE 
vs IRAS, and the solid line is the average of the two. 
Computing both linear regressions allows us to conclude 
that the trend observed here is not due to noise. The 
effect of detector noise is clearly seen at small scale and at 
low brightness where the dispersion of the slope values is 
increasing. There is also an increase of the dispersion at 
large scales (very clear at 60 /im), due to the large scale 
drift of the detector responsivity or to a bad zodiacal 
light correction. These two effects produce an additive 
effect with spatial variations at scales larger than a few 

^ Note here that, as stated at the end of this section, < > is 
the median value of the DIRBE data to limit the effect of variations 
of the IRAS zero level 
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Fig. 4. — Illustration of the use of the wavelet decomposi- 
tion to compute a gain difference between two maps at different 
scales. Here the mi and m2 are identical fBm images (shown in 
the top row) with a constant multiplicative factor between the two 
(rri2 = 1.389 X mi). The four bottom rows shows the wavelet 
decomposition of mi and m2. The right column shows the corre- 
lation between m2 and mi for the whole map and at each scale. 
The slope a of the linear regression is given in each panel. 



degrees, affecting predominantly the slope measured at 
larger scales and at low brightness. 

The curves in figures and [7| can be resonably well 
fitted by the following fonctional: 

Ro + A< h> +B < Ix>^ . (3) 

The values of _Ro: ^ and B, which all depend on A and 
scale I, are given in Table [3 The brightness < > is 
the median value of the DIRBE data, to limit the effect 
of variations of the IRAS zero level. 

4.4. Reponsivity correction 

In the previous sections we showed that the IRAS re- 
sponsivity to diffuse emission varies with scale and, at 60 
and 100 /im, also with brightness. At these wavelengths 
we were able to establish an analytical form of the re- 
sponsivity variations (Eq.|31) by comparing the IRAS and 
DIRBE data. This expression is also valid for the 12 and 



Fig. 5. — Example of a wavelet decomposition of a typical 60 ^im 
image. See Fig.|3]for details. 



25 ^im bands; in this case it reduces to R = Rq as the 
responsivity does not depend on the brightness {A = 
and B = 0). 

The responsivity correction to apply to the ISSA plates 
is the inverse of the expression given in Eq. I^KG = 1 /R) . 
In Fig. |S1 we show the responsivity correction factor as 
function of scale for the four bands. For the 60 and 
100 ^m bands we present three curves to show the vari- 
ation of the responsivity correction with brightness. 

The implementation of such a correction is made diffi- 
cult by the fact that it depends on brightness (at 60 and 
100 /im). We need to apply a correction that depend on 
scale but also on position to take into account brightness 
variations across the map. To achieve this we use a local 
approch based on a wavelet decomposition of each ISSA 
HCON. 

To avoid discontinuity problems between adjacent 
plates each HCON is put in a 600 x 600 map and the 
undefined regions are filled using other HCONs or, when 
IRAS data are unavailable, by DIRBE (properly scaled 
to take into account the gain difference between the two 
data sets). For each HCON at each wavelength, we first 
separated point sources from diffuse emission (see § ^ . 
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Fig. 6. — Responsivity of the 60 fim IRAS data compared 
to DIRBE as a function of the 60 /im DIRBE brightness and for 
four different scales. The dotted hnes are the slope of the linear 
regression of IRAS vs DIRBE and DIRBE vs IRAS. The solid line 
is the average of the two dotted lines. The dashed line is the fit to 
the average slope using Eg. O 



Fig. 7. — Responsivity of the 100 /im IRAS data compared to 
DIRBE as a function of the 100 ^m DIRBE brightness and for four 
different scales. See Fig. |H] for details. 



To correct the diffuse emission map we decomposed it 
(using the "a trou" algorithm) for the first 6 scales (cor- 
responding to 3', 6', 12\ 24', 48', f.6°, 3.2°). Then we 
computed the proper gain correction factor at each scale 
and for each pixel of the map using Eq. |2| To compute 
this correction at 12 and 25 fj,m bands we simply use the 
responsivity values given in Table |21 and interpolate at 
the 6 corresponding scales (from 3' to 3.2°). At the 60 
and 100 ^m we do the same but we also use the DIRBE 
brightness to compute the brightness dependance part of 
the responsivity correction. 

The wavelet coefficients of the diffuse emission map 
were then recombined to form the responsivity cor- 
rected diffuse emission map. Finally the point sources 
map, properly calibrated by applying the corresponding 
DC/ AC factor (see Table ^ was added to the corrected 
diffuse emission map. The responsivity correction ap- 
plied here does not correct the emission at scales larger 
than 3.2° degrees. We take care of this in the next and 
step of the processing. 

4.5. Zero level correction 

As it was mentioned earlier the zero level of the IRAS 
data is not well determined, due to detector drifts and/or 
residual zodiacal emission. To correct this large scale 
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Fig. 8. — Responsivity correction factor (G = 1/i?) for the 
four bands as a function of angular scale. For the 60 and 100 iivn 
bands the gain factor is given for three values of mean brightness. 
The solid, dotted and dashed lines are respectively for 1, 10 and 
50 MJy/sr at 100 fim and 0.3, 3 and 15 MJy/sr at 60 ^tm. The 
shaded grey area represents the angular scales on which DIRBE 
data are used to scale the IRAS data. 
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Fig. 9. — Plate number 81 at 12 fj,m. Comparison of the ISSA 
(top-left) and IRIS (top-right) maps. The difference between the 
two (ISSA-IRIS) is shown in the bottom-left panel. The ratio of 
the power spectra of the IRIS and ISSA maps is shown in the 
bottom-right panel. 
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Fig. 10.— Plate number 81 at 25 fim. See Fig. [HI 



responsivity and the spatially varying zero level of the 
ISSA maps we again used the DIRBE data that were 
better calibrated and that used a better zodiacal emis- 
sion model. Here we make the assumption that the zero 
level does not vary on scales smaller than 1° which is 
reasonable as detector drifts and zodiacal emission are 
not expected to vary on such scales. 
As for the responsivity correction we computed the 
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Fig. 11.— Plate number 81 at 60 fim. See Fig.lHl 
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Fig. 12.— Plate number 81 at 100 /.tm. See Fig.|2l 



offset map for every HCON, enlarged to 600 x 600 to 
make sure that we keep the continuity between adjacent 
ISSA plates. Empty pixels were filled with the average 
map (HCONO) , previously corrected for the responsivity 
variation, or DIRBE data when necessary. The HCON 
map was then convolved by the DIRBE beam and by a 
Gaussian beam (FWHM=40') to smooth out the DIRBE 
pixelisation. The DIRBE data were projected on this en- 
larged HCON map and convolved by the same Gaussian 
beam. At that point both data sets are on the same 




O.D10 0,100 

Fig. 13.— Top: Typical 100 urn IRIS image (left) together with 
its noise map (right). These two maps are in MJy/sr. Bottom: 
Coverage map (left) showing the number of data points averaged 
at each position (Ntot{x,y))- The power spectra of the 100 fim 
map and of its noise map are shown in the right panel. 



grid and at a 1° resolution'^ . We then computed the 
difference between the two maps (DIRBE-IRAS). This 
difference map was added to the HCON. When all the 
3 HCONs were corrected we computed the average map 
(HCON-0) using the appropriate coverage maps. This is 
the final step of the processing. 

5. PROPERTIES OF THE IRIS PLATES 

The processing of the 4893 HCON plates was done 
using the McKenzie cluster of the Canadian Institute 
for Theoretical Astrophysics. This Beowulf cluster is 
composed of 512 2.4 GHz Xeon processors with a total 
of 256 GB of RAM. We effectively used only 40 CPUs of 
the cluster at the same time and the total computation 
time for the whole processing (deglitching, destriping, 
responsivity and offset correction) was ^--^20 hours. 

FiguresEltoEl present the ISSA and IRIS plate num- 
ber 81 at the four wavelengths. Are also shown the dif- 
ference between the ISSA and IRIS maps and the ratio of 
their power spectra. These figures show that the destrip- 
ing described in 5 13. 21 removed very efficiently the stripes 
at all wavelengths. The second more noticeable improve- 
ment is the modification of large scale structures in each 
maps, due to the responsivity and zero level corrections. 
These two corrections removed very well residual zodi- 
acal light and detector drifts and improved significantly 
the calibration of the data. The correction of the large 

^ It is worth mentioning that point sources were removed in 
the IRAS and DIRBE data to compute the offset. Even after the 
gain correction, which reconciliate the calibration of point sources 
and diffuse emission, we had to remove point sources for the offset 
computation as bright sources at 12 and 25 /im have relatively 
strong tails that increase significantly their flux when convolved 
by the DIRBE beam. 



scale structures of the maps (that is due both to respon- 
sivity and zero level corrections) is clearly seen in the 
ratio of the power spectra. 

One interesting property seen on Figures |51 and 1101 
is that, even after the destriping, the power spectra 
ratio of the 12 and 25 /im images is 1 on scales 
between 5' and 1°. This is due to the responsivity 
correction applied (see Fig. (H)) that is higher than 1.0 at 
those scales. The noise reduction resulting in improved 
signal-to-noise ratios obtained from the destriping is 
almost counterbalanced by the multiplication of all the 
small scale fluctuations (including noise) by a factor 
greater than 1. Nevertheless the noise properties of 
the IRIS plates are much more easy to handle as they 
are now homogeneous. On the other hand, at 60 and 
100 /im the responsivity correction is not as important 
and there is a noise reduction that is about 15-20% at 
scales 5-15'. 

Finally, notice that the areas of the sky not covered by 
the IRAS survey are left empty in the IRIS product. 

5.1. Noise 

Following the i nethod described by 

iMiville-Deschenes et all l|20n2bD one can estimate 
the power spectrum of the noise of each IRIS map. 
To do so we take advantage of the fact that each final 
IRIS map is composed of a combination of up to three 
HCONs. Here we argue that the difference between two 
HCONs is dominated by the instrumental noise and that 
the power spectrum of such a difference map is equal to 
the power spectrum of the noise, multiplied by a scaling 
factor that depends on the coverage maps. 

If we make the assumption that the noise level B of 
the IRAS detectors is stationary, and if we neglect the 
photon noise (i.e. the noise level is dominated by the 
readout noise), the noise level at position {x,y) of the i*'' 
HCON map is 



B 



(4) 



where Ni{x,y) is the coverage map of the i HCON 
map. An example of a typical coverage map is shown 
on Fig. 1131 

The noise level of the co-added map (HCON-0) is 



CTo(a;,y) 



B 



^/Ntot{x,y) 



(5) 



where Ntot{x,y) = Ni{x,y) +N2{x,y) + N3{x,y). Simi- 
larly the noise level of the difference between two HCON 
maps (i and j) is 



(6) 




N,{x,y) Nj{x,y)' 



Combining the last two equations, the noise level of the 
co-added map becomes: 



cro{x,y) 



N,{x,y)Nj{x,y) 



Ntot{x,y){Ni{x,y) + N,{x,y)) 



(7) 



11 



To estimate the power spectrum of the noise of a given 
IRIS plate we first compute the difference map between 
two HCONs. The difference map is then normahzed fol- 
lowing Eq. [7| to take into account the spatial coverage 
of every HCON. An example of such a normalized dif- 
ference map is shown in Fig. 1131 together with its power 
spectrum and the power spectrum of the corresponding 
IRIS plate. At small scale (high k) these two power spec- 
tra meet perfectly which validates our processing. 

It is important to point out that the difference map 
used here to compute the power spectrum of the noise is 
by no mean a noise level map. It does not provide a way 
to estimate the noise level of every pixel in the map but 
it is only used to estimate the global properties of the 
noise (i.e. power spectrum, average level). 

The average noise level of an IRIS map can be esti- 
mated by taking the standard deviation of the difference 
map. Typical noise level values are given in Tabled for 
the four bands. Within 10% the IRIS noise level is the 
same as the ISSA noise level. This is due to the fact 
that the power spectrum of the noise is not flat (its slope 
varies from -1.5 to 0). Therefore the standard deviation 
of the difference map may be dominated by large scale 
structures in the noise. In such case, integral of the ISSA 
and IRIS noise power spectra on selected scales allows a 
better comparison. For large scales (1-10°), the noise 
level in the IRIS plates is about 30% smaller than that 
of the ISSA plates, for the four wavelengths. At smaller 
scales, as explained in fJJl we gain in sensitivity only at 60 
and 100 /im. The gain is about 15-20% for scales 5-15'. 

5.2. Effective angular resolution 

Due to the geometry of the IRAS detectors, the IRAS 
data have a higher spatial resolution in t he in-scan direc- 
tion than in the cross-scan (see TableP). iWheelock et~aLl 
|l993) have computed the PSF on few (4 to 5) point 
sources in each band. They showed that the PSF is non- 
circular due to difference in the in-scan and cross-scan 
resolution. They conclude that the angular resolution of 
the ISSA plates is between 4 and 5 arcmin in each band. 

Following Mivillc-Dcschcncs et al. ( 20021j) we used the 
power spectrum of ISSA plates to estimate the effective 
Point Spread Function (PSF) of the ISSA. If we make 
the assumption that the effective PSF is described by a 
Gaussian, the power spectrum of an image dominated by 
interstellar fluctuations can be expressed by the following 
equation: 

P{k) = exp (^^^ X Ak^ (8) 

where Ak^ describes the cirrus emission power spectrum 
(jGautier et al.l 11992^. We have done such a fitting for 
several IRIS plates at the 4 wavelengths. We have se- 
lected plates with relatively bright interstellar emission 
(not to be affected by noise) from which we have removed 
the bright point sources. The effective resolution of the 
four IRIS bands are given in Tabled We have checked 
that these values are compatible with Gaussian fitting on 
point sources. 

5.3. Calibration uncertainties 

At 60 and 100 /xm, the responsivity correction we ap- 
plied at scales smaller than 1.25° depends on an inter- 
polation between the DC/ AC factor (responsivity at 5' 



scale) and the responsivity at scale of 1.25° (deduced 
from the comparison of IRAS and DIRBE with scale and 
brightness - see S I4.3|l . The largest responsivity differ- 
ence between these two scales is for bright regions where 
it reaches ^20% (see Fig.lHJ. We estimate that the typi- 
cal uncertainty on the responsivity at scales smaller than 
1.25° is ~5%, for the four bands. 

At scales larger than 1.25° the IRAS data are forced 
to match the DIRBE data by applying an offset map 
(see i) I4.5|l . The application of such an offset map takes 
care of responsivity variations at scales larger than 1.25° 
but it also sets the zero level correctly. The uncertainty 
on the responsivity at large scales is thus limited by 
the DIRBE gain uncertainties which varies from 5.1 to 
15.1% (see Table nj. Considering the relatively small 
uncertainty of the responsivity at sub-degree scales 
we estimate that the overall uncertainty on the IRIS 
calibration is dominated by the uncertainty on the 
DIRBE data (including the uncertainty on the zodiacal 
model). 

In the comparison of the IRAS and DIRBE data, an 
additional uncertainty comes from the slight differences 
of the system spectral responses. To take into account 
these differences, the DIRBE data should have been mul- 
tiplied by the ratio of the IRAS to DIRBE color correc- 
tion prior to any comparison with the IRAS data (see 
§ Ib)). This ratio depends on the details of the wave- 
length response of the system but also on the shape of 
the intrinsic energy distribution that is unkown for each 
sky direction (especially at short wavelengths). No cor- 
rection was therefore applied to the DIRBE data. This 
induces an additional uncertainty at 100, 60, 25 and 12 
microns that is less than typically 6, 10, 3 and 10% for 
the diffuse medium respectively. 

5.4. Comparison with Schlegel et al. 

A first attempt to pro duce well calibrated 100 fj,m 
IRAS maps was done bv iSchlegel et all 1)1998(1 . These 
authors improved the destriping and zero level of the 
ISSA plates in a similar fashion as we did. One of the 
gr eat improvement of t he IRIS 100 /im data compared 
to ISchlegel et all ^98) is that we take into account the 
variation of the responsivity with brightness at scales 
smaller than 1° (while they apply a con stant factor of 
0.87). In addition the product given by ISchlegel et alJ 
l( 19981) does not include point sources. 

In Fig. 1141 we present a comparison of a 100 /im 
250 X 25 pixels sub- im a ge of the same region in ISSA, 
IRIS and ISchlegel et alJ l|1998D . One striking feature of 
this figure is the striping that is c ompletely removed in 
the IRIS and ISchlegel etaLl llT99l images. One feature 
of the IRIS product is also the fact that it includes point 
sources. One would also note the difference in angular 
resolution between the IRIS and pchlegel et al.l l)1998fl 
data. Contrary to ISchlegel et al.l l| 1998(1 we have taken 
great care to keep the full ISSA resolution. The difference 
in angular resolution between IRIS and ISchlegel et alJ 
(^98) is well seen in Fig.^|where we compare the power 
spectrum of a typical 100 /im map^. There is a very good 
agreement between these two data sets at scales larger 

Point sources have been subtracted from the IRIS map prior 
to compute the power spectrum. 
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Fig. 14. — Comparison of a ISSA, IRIS and Schlegel ct al. 250 X 250 sub-image at 100 fim. All the maps are in MJy/sr. The IRIS map 
combined the higher angular resolution of the ISSA map and the destriping and zero level correction of the Schlegel et al. map. In addition 
the IRIS map was corrected for the variation of the IRAS responsivity with scale and brightness and it also includes properly calibrated 

Compared to th e pre vious products (namely the 
ISSA and Schlege l "etall (0.998)), we made a major 
improvement of the calibration and zero level at all 
wavelengths (but especially at 12 and 25 /im) by com- 
paring with DIRBE and remove significantly the stripes 
using an eflicient iterative process that involve a wavelet 
decomposition of the map. Considering the relatively 
small uncertainty of the IRAS responsivity at sub-degree 
scales we estimate that the overall uncertainty on the 
IRIS calibration is dominated by the uncertainty on the 
DIRBE data (including the uncertainty on the zodiacal 
model). One other great interest of the IRIS maps is 
that they contain all the IRAS small scale information 
(including point-like sources) properly calibrated. With 
IRIS, we have reached the limit of the post-processing of 
the ISSA plates. Improvements over the IRIS product 
will be only possible by working directly on the raw 
IRAS data. 
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Fig. 15. — Comparison of the power spectrum of a typical 
100 fim image with IRIS and Schlegel et al. Left: Power spectra 
of the IRIS (point sources subtracted) and Schlegel et al. plates. 
Right: Ratio of the IRIS and Schlegel et al power spectra. The 
two plates have a very similar power spectrum except at high-A; 
due to the higher angular resolution of the IRIS plate. 



than ^ 30' but there is a significant difference at smaller 
scales due to the difference in the responsivity correc- 
tion a nd also to the lower resolution of the Schlegel ct al. 
lll99l data (6.1'). 

6. CONCLUSION 

We have built a new generation of IRAS maps that are 
rescaled to an absolute calibration coherent with both the 
DIRBE and the IRAS Point Source Catalog photometry. 
These maps are also much less contaminated by stripes 
than the previous IRAS images. This new generation of 
IRAS maps is named IRIS for Improved Reprocessing 
of the IRAS Survey. The IRIS data are available at 
Ihttp:/ /www.cita.utoronto.ca^mamd/IRIS and 
[http:// www. ias.u-psud.fr /iris They can also be 
obtained from http://www.ipac.caltech.edu and 
through the Centre de Donnees astronomiqucs de Stras- 
bourg (CDS - http://www.cdsweb.u-strasbg.fr) via 
Aladin. The data products are detailed in Appendix lO 



The IRIS data will be of great use for the whole com- 
munity. In the context of the numerous ongoing and 
future infrared/submm/mm missions (Spitzer, Herschel, 
Planck, JWST, ASTRO-F...), they will be very helpful to 
prepare and analyse the observations. In particular, with 
an angular resolution comparable to the high frequency 
bands of Planck, these data will be very complementary, 
providing unique dust templates. These data can also be 
used by themselves to conduct still original science. As 
an example, in a forthcoming paper, we will present a 
study of the properties of the cirrus dust emission. 
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APPENDIX 
POINT SOURCE IDENTIFICATION 

One important part of the post-processing of the ISSA maps presented here is related to the identification and 
removal of point sources, or more generally of small scale structures. The removal of point sources is used at several 
steps in the processing. As the calibration of the diffuse emission and of point sources are not the same in the ISSA 
plates, they have to be treated separately. Point sources identification is also used for the deglitching and destriping 
of the HCONs and for the determination of the zero level map. 

The difficulty here is to develop a point source identification method that is equally efficient in very diffuse regions 
where the noise is significant and in crowded Galactic plane fields. The poor-man method to identify point sources is to 
apply a boxcar filtering to the image and to look for highly deviant small scale fiuctuations. For images with relatively 
homogeneous background the histogram of small scale brightness fluct uation is dominated by a quasi-Gaussian core, 
partly due to noise, with a positive wing due to point sources (see Fig. IA16|) . In this case point sources can easily be 
identified as beeing all the pixels with fluctuation values higher than na where a is the width of the Gaussian core 
and n is a given threshold. 

For images with strongly varying diffuse emission this method does not work as well. In Fig. IA16I (middle row) we 
present the identification of point sources with this traditional sigma-clipping method. In this case the shape of the 
small-scale fluctuation histogram is very different and a much larger fraction of the pixels in this image have a small 
scale fluctuations above na, where a is the width of the central component of the histogram. The result (middle-right 
panel) is that most of the structure in the bright parts of the image has been identified as point sources. 

For bright regions the small-scale fluctuations are not dominated by noi se fluctuations . The core of the small-scale 
fluctuation histogram comes from the diffuse emission itself. As shown bv lGautier et al.l ljl992j) . the amplitude of the 
diffuse emission fluctuations P{1) at a given scale / increases with the average brightness Bq: 

P{1) (X (Al) 

iGautier et al.l l)1992D gave a — 1.5. This indicates clearly that the fluctuations level of the diffuse emission increases 
with brightness. This is clearly seen in Fig. IA161 where bright regions have generally stronger small scale fluctuations. 

To take this effect into account we rather use the contrast of the small scale fluctuations (i.e. fluctuations / brightness) 
to identify point sources, instead of the small scale fluctuation level itself. This way we take into account the variation 
of the diffuse emission fluctuations with brightness (see Eq. IA1|) . The result of this method is shown on the bottom 
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Fig. A16. — Exemple of point source identification. 



row of Fig. IA16I much less pixels were clipped. This method is much more flexible than the traditional sigma-clipping 
method as it adapts itself to the background level. 

At this point it is important to mention that the contrast method can't be applied in diffuse region where noise 
contributes significantly to the small scale brightness fluctuations. For each band we used the traditional sigma-clipping 
for regions with boxcar filtered brightness lower than 3, 3, 4 and 4 MJy/sr at 12, 25, 60 and 100 /xm respectively. For 
regions with brightness higher than those values we used the contrast method to identify point sources. 

IRAS AND DIRBE FILTERS 

Following the IRAS convention, the DIRBE spectral intensities l^, are expressed in MJy/sr at fixed nominal fre- 
quencies assuming the source spectrum is ^'Ij,=constant (i.e. constant intensity per logarithmic frequency interval). 
Since the source spectrum has not a constant intensity per logarithmic frequency interval, a color correction has to be 
applied to obtain an accurate intensity. The color correction factor is defined such that; 

I^g{actual) — I^g{quoted)/cc (Bl) 

where 1^,^ (actual) is the actual specific intensity of the sky at frequency vq, Ii,g{quoted) is the corresponding value 
given with the IRAS or DIRBE convention and vq is the frequency corresponding to the nominal wavelength of the 



15 



Wavelength 


(T (K), /3) 




100 fim 


('14 2~) 


1.1019 




flfi 2~) 


1.0737 




(18 2) 


1.0525 




(20, 2) 


1.0354 






1.0211 




(14, 1.5) 


1.1146 




(16, 1.5) 


1.0852 




(22, 1.5) 


1.0321 




{-, -) 


1.0510 


60 /^m 


(40 1) 


1.0476 




(60, 1) 


0.9452 




(40, 0) 


1.1187 




(60, 0) 


0.9947 




(-, -) 


1.1045 


25 


(- 0) 


0.8941 




(-. 1) 


0.8036 




(-, 2) 


0.7265 




(-, -) 


0.9728 


12 /im 


(-■ 0) 


1.0868 




(-■ 1) 


1.1723 




(-, 2) 


1.2552 




(-, -) 


1.0938 



TABLE B3 

Ratio of IRAS to DIRBE color corrections for spectra in the form Ii, oc !^''B^(T). When no values are given for the 

TEMPERATURE, THE CONSIDERED SPECTRUM IS IN THE FORM 1^ (X 1^^ . WHEN NO VALUES ARE GIVEN FOR BOTH THE TEMPER ATURE AND 
SPECTRAL INDEX, THE CONSIDERED SPECTRUM IS THE STANDARD DIFFUSE ISM SPECTRUM FOR 10^" At/CM^. 



considered band. With these definitions: 

^ J (Iu/Iuq) actual Rydv ^^^-j 

/ {va/v)Ri,dv 

where {lul Ivo)actuai is the actual specific intensity of the sky normah sed t o the intensity at frequency vq and is 
the system spectral response. The spectral responses are given in Fig. IB 171 for the four IRAS and DIRBE bands. 

When we compare the DIRBE and IRAS intensities, we compare Iua{quoted)DiRBE and I^^-,{quoted)iiiAS- For 
the same point on the sky and in the same beam (i.e. when IRAS data are convolved by the DIRBE beam), 
I (actual) DIRBE = Ii^a (actual) I HAS- This leads to: 

(quoted) IRAS cciras ,-oo\ 

- (a6) 



I uo(quoted) dib.be ccdirbe 

Therefore when comparing the DIRBE and IRAS data to calibrate the IRAS data, we should correct the DIRBE data by 
the ratio of color corrections. This ratio is varying across the sky. It is quite easy to compute it in the Far-Infrared at 
large scales where we can infer the dust emissivity and temperature using DIRBE data (using an iterative procedure). 
At smaller scales and at 12 and 25 /.tm, it is impossible to compute this ratio since we do not know the exact shape of 
the intrinsic energy distribution for most of the sky. As an example, ratios of color corrections are given for different 
spectral energy distribution in Table lB3l 

DATA PRODUCTS 

The "standard" IRIS data ("standard" means for the processing described in this paper) are available at 
|http: / /www.cita.utoronto.ca'^mamd/IRIS and http: / /www. ias.u-psud.fr/iris On these two sites the data are available 
in the original ISSA format: 430 maps of 500 x 500 pixels for each band. I'he individual HCON can also be downloaded 
from these sites. In the near future we expect to provide Healpix version of these maps on these two sites. We also 
expect to provide versions of the maps without point sources and with an improved zodiacal light subtraction. 
The IRIS data product is also available through the Centre de Donnees astronomiques de Strasbourg (CDS - 



http : / / www . cdsweb . u-strasbg . fr ) via Aladin and at the Infrared Processing and Analysis Center (IPAC). 
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Fig. B17. — System spectral responses for the IRAS (continuous lines) and DIRBE (dashed lines) 100, 60, 25 and 12 /im bands (Figs, 
a. to d. respectively). 



